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ABSTRACT 

In  automobile  frontal  impact,  for  given  restraint 
characteristics  and  prescribed  impact  speed  and  crash 
deformation,  what  is  the  optimal  vehicle  crash  pulse  that 
produces  the  lowest  peak  occupant  deceleration?  In  this 
paper,  based  on  a  lumped-parameter  model  of  the 
occupant-vehicle  system,  the  problem  is  treated  as  a 
best  disturbance  problem  of  the  optimal  protection  from 
impact.  The  optimum  kinematics  of  the  occupant  in 
frontal  impact  is  found.  For  linear  restraint 
characteristics,  the  theoretical  optimal  crash  pulse  is 
obtained,  and  the  relation  of  the  peak  occupant 
deceleration  to  the  impact  speed,  crash  deformation, 
and  vehicle  interior  rattlespace  is  established.  The 
optimal  crash  pulses  for  passive,  active,  and  pre-acting 
restraint  systems  are  discussed.  Numerical  optimization 
is  formulated  to  find  the  optimal  crash  pulse  for  general 
restraint  systems.  The  occupant  responses  under  the 
constant-deceleration  crash  pulse  and  the  half-sine 
crash  pulse  are  analyzed  for  linear  restraint  systems. 
Comparisons  are  made  between  the  optimal  crash  pulse 
and  these  two  “non-optimal”  crash  pulses. 

INTRODUCTION 

The  structural  crashworthiness  plays  important  roles  in 
the  prevention  and  reduction  of  occupant  injuries  in 
automobile  crashes.  To  attain  the  optimal 
crashworthiness  of  an  automobile  structure,  optimization 
is  an  effective  tool  and  becomes  a  common  practice  in 
the  vehicle  design.  Optimization  can  now  be  conducted 
in  a  virtual  environment  based  on  computational 
modeling,  simulation,  and  analysis  of  a  system. 

The  optimization  problem  of  the  vehicle  crashworthiness 
can  be  dealt  with  analytically  with  two  approaches,  direct 
and  indirect  [1],  The  direct  approach  uses  physical 
parameters  of  the  vehicle  structure  as  the  design 
variables  in  optimization.  The  indirect  approach  handles 
the  issue  by  a  two-step  strategy:  first,  for  the 
minimization  of  the  occupant  injuries,  find  the  optimal 
response  for  the  vehicle  structure;  then,  solve  an 
inverse  design  problem  for  the  vehicle  structure 
according  to  its  optimal  response  or  use  the  optimal 
response  as  a  guideline  for  the  vehicle  structure  design. 


The  study  in  this  paper  will  take  the  indirect  approach 
and  deal  only  with  the  problem  in  the  first  step  for  frontal 
impact,  i.e.,  the  optimal  vehicle  crash  response.  In 
particular,  the  issue  to  be  addressed  is:  in  frontal  impact, 
for  given  restraint  characteristics  and  prescribed  impact 
speed  and  crash  deformation  of  the  vehicle,  what  is  the 
optimal  vehicle  crash  pulse  for  mitigating  occupant 
injuries? 

A  number  of  studies  have  been  performed  on  the 
vehicle  crash  pulse  in  the  last  three  decades  or  so. 
Several  late  studies  were  particularly  concerned  with  the 
problem  of  the  optimal  crash  pulse  [1-4],  In  [2],  based 
on  a  one-mass  model  and  a  rigid  multi-body  model,  the 
vehicle  crash  deceleration  curve  was  divided  into  finite 
segments  and  their  effects  on  the  occupant's  injury  were 
investigated  quantitatively  using  a  sensitivity  analysis. 
The  optimal  deceleration  curve  for  minimizing  the 
occupant’s  injury  was  obtained  by  iterative  calculations. 
In  [3],  based  on  a  two-mass  one-dimensional  model, 
pure  rectilinear  form  is  assumed  as  the  ideal  pulse  for 
the  occupant  deceleration.  Three  functions  were  formed 
to  approximate  the  rising  edge  of  the  rectilinear  form 
and  corresponding  initial  parts  of  the  vehicle  crash  pulse 
were  derived.  In  [1]  based  on  a  spring-mass  model,  the 
problem  of  the  optimal  crash  pulse  is  addressed  by 
applying  an  energy  relationship  between  the  vehicle 
crash  pulse  and  the  occupant  deceleration  in  a  domain 
of  the  relative  motion  between  the  vehicle  and  the 
occupant.  It  was  found  that  for  prescribed  impact 
velocity  and  vehicle  crash  deformation,  the  optimal 
pulse  is  one  that  consists  of  an  impulse,  a  subsequent 
zero-acceleration  period,  and  finally  a  constant  level 
period.  In  [4],  the  problem  was  also  described  by  a  one- 
mass  model,  the  crash  pulse  was  discretized  with 
respect  to  the  crash  deformation,  and  a  solution 
procedure  using  numerical  optimization  was  proposed  to 
find  the  optimal  crash  pulse. 

In  this  paper,  the  problem  of  the  optimal  vehicle  crash 
pulse  is  treated  as  a  best  disturbance  problem  of 
optimal  protection  from  impact. 

SYSTEM  MODELING  AND  PROBLEM 
STATEMENT 


The  models  used  for  computational  modeling  of 
automobile  crashworthiness  can  be  classified  as 
lumped-parameter  models,  rigid  multi-body  (RMB) 
models,  finite  element  (FE)  models,  and  integrated 
models  [5],  These  models  provide  different  levels  of 
abstraction  of  the  problem,  and  thus  have  different 
applications.  In  frontal  impact,  the  occupant  can  be 
subject  to  injuries  to  various  regions  of  the  body,  such 
as  the  head,  thorax,  upper  extremities,  and  lower 
extremities.  The  injury  criteria  for  measuring  these 
injuries  are  usually  expressed  by  impact  responses 
(such  as  accelerations)  in  respective  regions.  When  the 
prevention  and  reduction  of  injuries  to  the  occupant  are 
considered  in  general,  the  occupant  can  be  reasonably 
treated  as  a  point  mass  and  its  peak  acceleration 
(deceleration)  can  be  used  as  the  injury  criterion.  When 
the  vehicle  crash  pulse  is  under  investigation,  the 
vehicle  is  treated  as  a  whole  body  and  only  its  gross 
motion  is  considered. 

Therefore,  automobile  frontal  impact  is  described  by  a 
lumped-parameter  model  shown  in  Fig.  1,  where  the 
occupant  is  represented  by  a  point  mass  m , 
x(t.)  describes  the  gross  motion  of  the  vehicle,  andj/(/) 

is  the  motion  of  the  occupant  relative  to  the  vehicle.  The 
interaction  between  the  occupant  and  restraint  systems?, 
such  as  seat  belt  and  airbag,  is  simply  represented  by  a 
spring  and  a  damper,  but,  in  general,  it  can  be  described 
by  a  control  force  u(y,y,t)  .  For  the  vehicle  motion, 
a(t)  =  -x(t),  (1) 

usually  referred  to  as  the  vehicle  crash  pulse,  and 
Dv  -  max  {*(7)} ,  (2) 

defined  as  the  vehicle  crash  deformation.  The  free 

space  between  the  occupant  and  the  vehicle  interior  is 
called  the  rattlespace  S0,  which  is  the  maximum 

allowable  space  for  the  occupant  to  move  forward  with 
respect  to  the  vehicle,  i.e., 

max{X0}^*V  (3) 

Since  this  space  has  important  effects  on  the  occupant 
response  also  and  is  a  major  factor  in  the  restraint 
system  design,  it  will  be  considered  in  addition  to  the 
crash  deformation. 


with  the  initial  conditions 

x(0)  =  0,  x(0)  =  v0 ,  (5) 

and 

_K°)  =  0,  y(0)  =  0  .  (6) 

where  v0  is  the  impact  velocity. 

The  question  to  be  addressed  in  this  paper  can  be 
expressed  as: 

Problem-1:  For  given  restraint  characteristics  and 
prescribed  impact  velocity  and  crash  deformation,  find 
the  optimal  crash  pulse  such  that  the  peak  occupant 
deceleration  is  minimized.  The  problem  can  be 
formulated  as:  Find  the  optimal  vehicle  crash 

pulse aQ(t) ,  such  that 

7,  (a0 )  =  min  {./,  (a)  \  J 2  (a)  <DV,J3  ( a )  <  S0  }  ,  (7) 

a 

where 

J,  =  max  {-[*(0  + KOI} .  (8) 

the  peak  occupant  deceleration, 

y2  =max{x(/)},  (9) 

the  maximum  vehicle  displacement,  and 

y3  =  max{X0).  (10) 

the  maximum  forward  excursion  of  the  occupant  in  the 
vehicle. 

From  the  standpoint  of  the  kinematics  of  the  occupant, 
Problem-1  can  be  reduced  to  a  more  generic  problem: 

Problem-2:  For  prescribed  impact  velocity,  crash 
deformation,  and  rattlespace,  find  the  optimal 
kinematics  of  the  occupant  such  that  the  peak  occupant 
deceleration  is  minimized  while  the  maximum  occupant 
forward  displacement  is  bounded.  The  problem  can  be 
formulated  as:  Find  the  optimal  kinematics  of  the 
occupant  w0(t)  ,  such  that 

yi(H'0)  =  min{y1(w)|y4(w)<£)o},  (ii) 

where 

y4=max{x(0+X0).  (12) 

t 

which  is  the  maximum  forward  displacement  of  the 
occupant,  and 

D0  =  Dv  +S0,  (13) 

the  allowable  maximum  forward  excursion  for  the 
occupant  in  an  inertial  frame.  Note  that,  unlike  a(t )  in 
Problem-1,  w(t )  is  not  a  control  function  of  the  system. 
Instead,  it  is  controlled  or  produced  by  a(l)  or 

*(y,y,  0- 

OPTIMAL  OCCUPANT  KINEMATICS 


The  equation  of  motion  of  the  system  is 

m(x+y)+u(y,y,t)  =  0, 


(4) 


2 


To  find  the  optimal  kinematics  of  the  occupant  for 
Problem-2,  according  to  duality  or  reciprocity  of 


optimization  [6],  the  dual  or  reciprocal  problem  of 
Problem-2  can  be  formulated  as  follows. 

Problem-3:  For  prescribed  impact  velocity,  vehicle 
crash  deformation,  and  rattlespace,  find  the  optimal 
kinematics  of  the  occupant  such  that  the  maximum 
occupant  displacement  is  minimized  while  the  peak 
occupant  deceleration  is  bounded.  The  problem  can  be 
formulated  as:  Find  the  optimal  kinematics  of  the 

occupant  w0(l) ,  such  that 

J 4 (wo )  ~  m111  (A  (w)  I  J\ (H0  -  Am) ,  (14) 

w 

where  Am  is  the  upper  bound  on  the  peak  occupant 
deceleration. 

To  find  the  theoretical  solution  of  Problem-3,  denote 
z(t)  =  x(t)  +  y(l),  (15) 

which  is  the  absolute  motion  of  the  occupant  with 
respect  to  an  inertial  frame  during  impact.  Consider  the 
motion  of  the  system  from  the  initial  of  impact  i  -  0  to 
the  instant  t  =  T0  when  the  occupant  comes  to  a  rest  for 


the  first  time.  That  is, 

v0  +  Jo  z(z)dr  =  0.  (1 6) 

The  velocity  of  the  occupant  z(t)  starts  with  v0  at  t  =  0 
and  decreases  to  0  at  /  =  T0 ,  i.e. , 
z(t)>0,  0<t<T0.  (17) 

Therefore,  forO</<F0,  the  displacement  of  the 
occupant  z(t )  ,  which  is  given  by 

z (t)  =  z(z)dz ,  (1 8) 

increases  monotonically  with  respect  to  time  and 
reaches  its  maximum  value  at  t  =  T0  ,  that  is, 

ma x{z(t)}  =  z(TQ)  .  (19) 

To  minimize  the  maximum  occupant  forward 
displacement  J4, 

J ,  =  max(x(f)  +  j/(f)}  =  max  {z(t)}  =  z(T0 )  ,  (20) 

t  f€[o.r0] 


T0  should  be  minimized.  According  to  the  definition  of 
T0  (Eq.  (16)),  if  the  deceleration  of  the  occupant  -  z(t) 
takes  the  maximum  allowable  value  A  ,  that  is, 

m=-Am,  (2d 

T0  is  minimized,  and 

T0  =  (22) 

Am 

This  means  that  in  order  to  minimize  the  peak 
displacement  of  the  occupant,  the  deceleration  of  the 
occupant  should  remain  constant  at  the  value  ofAm. 


This  is  the  optimal  kinematics  of  the  occupant  for 
Problem-3. 

Based  on  the  duality  or  reciprocity  between  Problem-2 
and  Problem-3,  the  optimal  kinematics  of  the  occupant 
for  Problem-2  can  be  stated  as:  In  order  to  minimize  the 
peak  deceleration  of  the  occupant,  the  deceleration  of 
the  occupant  should  remain  constant  at  a  value,  which 
is  denoted  as  Am  and  is  given  by 

Am=v2Q/(2D0).  (23) 

This  is  also  the  optimal  kinematics  of  the  occupant  for 
Problem-1. 

THEORETICAL  SOLUTION  OF  OPTIMAL  CRASH 
PULSE 

Linear  Restraint  Characteristics 
Optimal  Crash  Pulse 

Suppose  the  onset  of  the  vehicle  impact  is  at  t  =  0 . 
Consider  the  motion  of  the  occupant  right  after  the  onset 
of  impact  (/  =  0+ )  until  a  stop  for  the  first  time  (/  =  T0). 
As  the  initial  conditions  of  the  occupant  are 
z(0+)  =  0,  z(0+)  =  v0,  (24) 

the  optimal  kinematics  of  the  occupant  is 
2(0  = 

z(t)  =  v0-Aj,  0  ~<t<T0,  (25) 

2(0  =  v0t-±Amt2 

where  Am  is  yet  to  be  determined.  From  Eq.  (4) 

u(y,yJ)=mAm.  (26) 

Suppose  the  characteristics  of  a  restraint  system  can  be 
expressed  by  a  linear  spring  and  a  linear  damper.  That 
is, 


u(y,y,t)  =  ky+cy . 

(27) 

Then, 

ky  +  cy-  mAm , 

(28) 

with  initial  conditions 

^(o+)  =  o,ko+)  =  o. 

(29) 

The  solution  of  Eq.  (28)  is 

/  x  mA„,  ,, 

y(0  =  — — (l-e  e ) 

k 


mA„  — <  , 

y{t)  = - -e  c  ,  0+  <  t  <  70 .  (30) 

c 


From  Eqs.  (15)  and  (25),  the  optimal  motion  of  the 
vehicle  can  be  expressed  as 


3 


(40) 


km 

x(t)  =  -Am(\-  — e  (31  a) 

c 

m  -k-‘ 

x(0  =  v0  -Am(l+— e  c  ),  0  +<t<T0,  (31b) 

c 

*(0  =  V~i4/ - -^(l-ec).  (31c) 

k 

It  follows  that 

x(0+)  =  0,  x(0+)  =  v0  -  Amml  c  .  (32) 

Compare  to  Eq.  (5),  from  /  =  0  to  /  =  0+,  the  velocity 
change  of  the  vehicle  is 

Av  =  x(0)-x(0  +)  =  -Am.  (33) 

c 


To  produce  this  velocity  change,  the  acceleration  of  the 
vehicle  needs  to  have  an  impulse  I x  at  the  initial  of 

impact  {t-  0).  Theoretically,  this  impulse  can  be 
represented  by  a  Dirac  delta  function  8{t) ;  that  is, 

A  ffi 

J*(0  =  — *(0.  (34) 

c 

where  8{t)  is  given  by 


■5(0  =  j 


(35) 


and 


f  8{t)dt  =  1 .  (36) 


As  such,  from  the  onset  of  impact  {t  =  0)  to  the  first 
stop  (t  =  T0 ),  the  optimal  vehicle  acceleration  is 

...  \  .  r,  w  ,  km  -~i 

x(0  =  -AJl  +  -S(t) - c],0<f<ro.  (37) 

c  c 

Since  impulse  /r(f)  results  in  only  velocity  change  but 
no  displacement  change  of  the  vehicle,  Eq.  (31c)  still 
holds  for  0  <  /  <  T0 . 


Note  that  the  exponential  terms  in  Eqs.  (30),  (31),  and 
(37)  diminish  with  time.  After  a  sufficiently  long  time 
they  can  be  neglected. 


Determination  of  A 

m 

If  restraint  characteristics  are  such  that  the  maximum 
forward  excursion  of  the  occupant  relative  to  the  vehicle 
reaches  the  prescribed  rattlespace,  in  other  words,  if  the 
rattlespace  is  fully  used,  according  to  Eq.  (23), 

,.2 

(38) 


Am  2(DV  +S0) 

It  follows  from  Eq.  (31c)  that 

max{x(/)}  =  x(T0)  »  Dv  +S0 


mA . 


=  D„ 


where  the  exponential  term  is  neglected.  Then, 


where 


(41) 


m 


which  can  be  considered  as  the  natural  frequency  of  the 
restraint-occupant  system.  From  Eqs.  (38)  and  (40) 


<y2  = 


2S0(DV+S0) 


(42) 


which  represents  a  requirement  for  the  restraint 
characteristics.  This  means  that  when  the  natural 
frequency  of  the  restraint-occupant  system  satisfies  Eq. 
(42),  the  rattlespace  S0  will  be  fully  used  for  the 


occupant’s  excursion  with  respect  to  the  vehicle. 


In  general,  however,  restraint  systems  are  designed 
such  that 

max{v(0}<^0-  (43) 

i 

i.e.,  the  rattlespace  is  not  fully  used.  In  this  case,  it 
follows  from  Eqs.  (23)  and  (31c)  that 

2  j 

max{x(f)}  =  x(T0  )  *  - OLb. «  Dv ,  (44) 

'  2  A„  k 

m 

where  the  exponential  term  in  Eq.  (31c)  is  neglected 
also.  From  Eq.  (44), 


*Pv  (  L  ,  2mVl 

2m  ]j  kDj 


1) 


-1) 


(45) 


which  indicates  that  the  peak  occupant  deceleration  or 
the  amplitude  of  constant  deceleration  depends  upon 
the  impact  speed,  vehicle  crash  deformation,  and  the 
natural  frequency  of  the  restraint-occupant  system. 


According  to  Eqs.  (30)  and  (45),  Eq.  (43)  means  that 

CO  > - , 

2S0(D  +  S0) 


(46) 


which  represents  a  requirement  for  the  restraint 
characteristics  such  that  the  rattlespace  is  not  fully  used. 
It  can  be  inferred  that  if 


co  2  < 


2  S0(D  +  S0)' 


(47) 


the  excursion  of  the  occupant  will  exceed  the  space 
limit.  This  is  a  situation  not  allowed  to  happen. 


Illustration  of  Optimal  Crash  Pulse 
From  Eq.  (37),  the  optimal  crash  pulse  can  be 
expressed  as 


(39) 
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Ait ) 


.  M  m  km  - 

Am[\  +  -S(t ) - re 

c  c 


1 


1 


=  Am[\^-^dit)~—e  2‘ l  0<t<To 


2  £a> 


4C1 


where 


(48) 


<r=- 


„  ~-~7T=-  (49> 

2  com  2  4  km 

which  is  referred  to  as  the  damping  ratio  of  the  restraint- 

occupant  system.  Introducing  co  and  ^  reduces  the 

parameters  in  Eq.  (48)  and  avoids  providing  a  value  for 
m  ,  which  varies  from  occupant  to  occupant. 


For  the  illustration  of  the  optimal  crash  pulse,  choose 
the  values  for  the  parameters  in  Eqs.  (45)  and  (48)  to  be 
representative  of  automobile  frontal  impact  [1]: 

v0  =15.56  m/s  (35  mph),  Z)v=0.71  m  (28  in), 

co2  =42.43  rad/s,  and  ^  =  0.20.  The  optimal  crash 

pulse  together  with  the  velocity  and  displacement  of  the 
vehicle  and  the  optimal  kinematics  of  the  occupant  are 
illustrated  in  Fig.  2.  Note  that  the  impulse  of  the  delta 
function  at  t  =  0  should  lead  the  deceleration  of  the 
vehicle  at  /  =  0  to  go  infinite,  but  in  the  figure  only  a 
finite  large  value  can  be  displayed.  After  /  =  0,  the 
acceleration  of  the  vehicle  jumps  to  a  large  positive 
value  and  then  decays,  and,  after  a  certain  period  of 
time,  it  basically  remains  constant  at  a  negative  value. 
With  this  optimal  crash  pulse,  the  deceleration  of  the 
occupant  remains  constant  during  impact.  After  a 
certain  time,  the  relative  acceleration  and  velocity 
between  the  occupant  and  the  vehicle  become  zero,  and 
the  relative  displacement  between  them  remains 
constant. 


Trade-off  Curve —  A„  vs.  Sn 

m  0 

For  a  prescribed  optimal  crash  pulse  (v0  and  Dv  are 

given),  if  the  natural  frequency  of  a  restraint  system 
meets  the  requirement  of  Eq.  (42),  the  given  rattlespace 

S0  will  be  fully  used  and  then  the  peak  occupant 
deceleration  will  be  minimized.  For  the  restraint  system 
design,  S0  is  an  important  factor  and  can  be  used  as  a 

common  basis  for  the  comparison  of  various  designs. 
The  relationship  between  the  minimum  peak  occupant 

deceleration  Am  and  the  prescribed  rattlespace  S0  can 

be  expressed  in  the  form  of  a  so-called  trade-off  curve. 
According  to  Eq.  (38),  the  relationship  between  the  peak 
occupant  deceleration  Am  and  the  rattlespace  S0  is 

shown  in  Fig.  3  by  the  blue  curve  for  a  range  of  S0, 
where  v0  =15.56  and  Dv  =0.71.  The  curve  indicates 


that  as  the  rattlespace  S0  increases,  the  peak  occupant 


<  Time  (s) 


Time  (s) 

Figure  2.  The  motions  of  the  vehicle  and  the  occupant 


Figure  3.  The  relationship  between  the  peak  occupant 
deceleration  and  the  rattlespace 

General  Restraint  Characteristics 

In  general,  the  characteristics  of  restraint  systems  can 
be  so  complicated  that  they  cannot  be  reasonably 
described  by  a  linear  spring  and  a  linear  damper. 
However,  regardless  of  the  complexity  of  restraint 
characteristics,  once  they  are  given,  the  relative  motion 
of  the  occupant  with  respect  to  the  vehicle  can  be 
determined  by  solving  Eq.  (26),  and  then  the  optimal 
motion  including  the  optimal  crash  pulse  of  the  vehicle 
can  be  obtained.  Whereas  the  optimal  crash  pulse 
depends  upon  the  characteristics  of  particular  restraint 
systems,  certain  general  aspects  of  it  can  be 
ascertained  based  on  the  preceding  analyses  and 
results. 

•  If  restraint  systems  are  passive  and  have  certain 
damping,  an  impulse  is  required  at  the  initial  of  the 
optimal  crash  pulse. 
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•  If  restraint  systems  are  active  and  able  to  generate 
an  impulsive  action  on  the  occupant  actively,  an 
initial  impulse  may  not  be  required  for  the  optimal 
crash  pulse. 

•  The  major  portion  of  the  optimal  crash  pulse  after 
an  initial  period  tends  to  be  flattened. 

NUMERICAL  SOLUTION  OF  OPTIMAL  CRASH 
PULSE 

When  restraint  characteristics  are  nonlinear,  or  when 
restraint  characteristics  are  linear  but  the  motion  of  the 
vehicle  is  subject  to  other  constraints  such  as  a  limit  on 
the  amplitude  of  the  deceleration,  the  optimal  crash 
pulse  cannot  be  readily  obtained  analytically.  Instead, 
numerical  optimization  can  be  utilized  to  find  solutions. 

In  order  to  use  numerical  optimization  methods  to  solve 
Problem-1,  discretize  the  time  interval  [O,T0]  into 
identical  subintervals.  In  addition,  represent  the  crash 
pulse  a(t )  on  each  of  the  subintervals  by  a  constant 
value,  that  is, 

a{t)  =  aj  for  (l -i)At  <  t  <  i  At,  i  =  [50] 

where  N  is  the  number  of  subintervals,  At  is  the  length 
of  each  subinterval,  and  At  =  T0/ N .  Denote 

A  =  [a,a2...aN  f  [51] 

as  the  vector  of  design  variables.  The  motion  of  the 
occupant  relative  to  the  vehicle  is  discretized 
accordingly 

yi  =  y(i  At) ,  yi  =  y(i  At ) ,  i  =  1, . . . ,  N,  [52] 

and  can  be  obtained  from  the  numerical  integration  of 
Eq.  (4).  The  system  performance  criteria  of  Eqs.  (8)  - 
(10)  can  be  expressed  by  the  discretized  system 
responses: 


J\  (A)  =  max(-(.T 

l£j£A/ 

(53) 

J,  (A)  =  max{jc,}, 

1  \SiiN  ' 

(54) 

and 

■A  (A)  =  max  {>/,}. 

(55) 

Then  Problem-1  can  be  formulated  as  a  numerical 
optimization  problem: 

Design  Variables :  A  , 

Objective  Function'.  min{  J]  (A)  };  (56) 

Constraints.  (A)  <  D 

A  (A )<s0; 

al<a<au 


where  AL  and  Au  are  the  lower  and  upper  bounds  on 

the  crash  pulse.  A  nonlinear  optimization  method  can  be 
used  to  solve  this  problem  [7], 

OTHER  CRASH  PULSES 

The  optimal  crash  pulse  given  by  Eq.  (48)  is  meaningful 
only  in  a  theoretical  sense,  since  the  initial  impulse, 
which  is  represented  by£(/),  cannot  be  produced  in 

practice.  To  be  more  practical,  the  vehicle  deceleration 
during  impact  can  be  assumed  to  remain  constant;  that 
is,  the  crash  pulse  is  a  constant-deceleration  pulse. 
More  often,  the  vehicle  deceleration  in  frontal  impact  is 
represented  by  a  half-sine  pulse.  What  are  the  occupant 
responses  under  these  two  “non-optimal"  crash  pulses? 
The  issue  will  be  investigated  in  the  following  analyses. 
Again,  a  restraint  system  is  assumed  to  have  linear 
characteristics,  as  described  by  Eq.  (27). 

It  can  be  reasonably  assumed  that  the  peak  occupant 
deceleration  occurs  within  the  duration  of  the  vehicle 
crash  pulse.  Therefore,  the  occupant  responses  will  be 
considered  within  the  time  span  of  the  crash  pulse  in  the 
following  analyses.  In  the  following  computations, 
choose  impact  velocity  v0  =15.56  m/s  and  crash 

deformation  Dv  ~  0.71  m,  to  be  representative  of 
automobile  frontal  impact  [1]. 

Constant  Deceleration  Pulse 


Suppose  the  vehicle  deceleration  is  constant,  that  is, 


x(t)  =  -Av,0<t<Tv, 

(57) 

where 

A  -iL 

A'~2  D' 

(58) 

and 

II 

' 

(59) 

where  Tv  is  the  duration  of  the  crash  pulse.  The  velocity 


and  displacement  of  the  vehicle  are 

x(t)  =  v0-Avt,  (60) 

and 

*(0=V“4/2.  (61) 

From  Eq.  (28),  the  equation  of  motion  of  the  occupant 
with  respect  to  the  vehicle  is 

yF2Ccoy  +  a)2y~Av,  (62) 

with  initial  conditions 

MO)  =  0,  K0)  =  0 .  (63) 

By  solving  this  equation,  the  motion  of  the  occupant 
relative  to  the  vehicle  is  given  by 

y(.0  =  yc(t)+yp(t).  (64) 

Here, 


6 


>,(o= 4r>- e 


r  COS(7l“C2^  -«>-)]-  (65) 


where 


a  =  tan  '( 


<T 


0. 


and 


40  =  4  sm(^j]-C2 cot  +  fiy), 

where 


4  = 

and 


jCil  +Cly  +  2£CXyC2y 


Py  ~  tan' 
In  Eq.  (68), 
C  -—( 

ly  (O2 


i-C 


iy 


£C]y  +C2y 


Vi-? 


cos  av  -1) 


C2y=^f(may 


C 


CO 


a/i4 


cosa  ) 


(66) 


(67) 


(68) 


(69) 


(70) 


The  absolute  motion  of  the  occupant  is  the  superposition 
of  the  vehicle  motion  and  the  relative  motion: 
z{t)  =  x(t)  +  y(t) ,  (71) 

where  x{t)  is  given  by  Eq.  (61)  and  y(/)  by  Eq.  (64). 


The  occupant  responses  under  a  constant-deceleration 
pulse  are  illustrated  in  Fig.  4,  where  co  =  35.43  and 
It  is  shown  that  the  occupant  deceleration 
during  impact  is  not  constant,  which  means  that  the 
optimal  kinematics  of  the  occupant  is  not  attained.  In 
fact,  under  a  constant-deceleration  crash  pulse,  it  is 
impossible  for  a  passive  restraint  system  to  provide  the 
required  protection  to  the  occupant  such  that  the  optimal 
kinematics  of  the  occupant  is  achieved  [8],  However, 
the  restraint  characteristics  of  a  passive  system  can  still 
be  optimized  to  minimize  the  peak  occupant 
deceleration  as  much  as  possible.  Table  1  provides  the 
values  of  the  required  rattlespace,  optimal  frequency, 
and  minimal  peak  occupant  deceleration  corresponding 
to  each  level  of  damping  for  a  linear  passive  restraint 
system,  which  were  obtained  from  numerical 
optimization  [8], 

Based  on  the  results  given  in  Table  1,  the  relationship 
between  the  peak  occupant  deceleration  and  the 
rattlespace  is  displayed  in  Fig.  3  by  the  green  curve. 
The  curve  indicates  that  the  peak  occupant  deceleration 
increases  with  the  increase  of  the  rattlespace,  which 
seems  not  logical.  The  consequence  is  due  to  the 
damping.  In  the  range  of  [0.01  0.728],  larger  damping 
results  in  smaller  required  rattlespace  and  lower  peak 
occupant  deceleration.  However,  as  shown  in  Fig.  3,  for 
a  given  rattlespace,  the  peak  occupant  deceleration 


produced  by  a  constant-deceleration  pulse  is  much 
larger  than  that  resulting  from  the  optimal  crash  pulse. 
The  comparison  of  the  peak  occupant  decelerations 
between  the  optimal  crash  pulse  and  the  constant 
deceleration  pulse  is  shown  in  Table  2  for  a  series  of 
specific  rattlespaces. 


Time  (s) 


Figure  4.  The  motions  of  a  system  under  a  constant- 
deceleration  pulse 

Table  1.  Optimization  results  fora  constant  deceleration 
__ _ pulse  _ _ 


£ 

co  (rad/s) 

An  (»/»  ) 

Sn  (m) 

0.01 

34.63 

335.76 

0.280 

0.10 

34.83 

297.37 

0.243 

0.20 

35.43 

267.99 

0.207 

0.30 

36.47 

247.39 

0.176 

0.40 

38.10 

232.48 

0.147 

0.50 

40.60 

221.38 

0.120 

0.60 

44.61 

212.93 

0.094 

0.70 

52.73 

206.35 

0.064 

0.728 

60.39 

204.77 

0.048 

Table  2.  Comparison  of  peak  occupant  decelerations 
between  optimal  pulse  and  constant-deceleration  pulse 


So  (m) 

4  (m/s2) 

(Optimal 

Pulse) 

4  (m/s2) 

(Constant 

Deceleration 

) 

4.  -  4, 

4 

x  100% 

0.048 

159.71 

204.77 

28.2 

0.064 

156.40 

206.35 

31.9 

0.094 

150.57 

212.93 

41.4 

0.120 

145.85 

221.38 

51.8 

0.147 

141.26 

232.48 

64.6 

0.176 

136.63 

247.39 

81.1 

0.207 

132.01 

267.99 

103.0 

0.243 

127.03 

297.37 

134.1 

0.280 

122.28 

335.76 

174.6 

Half-sine  Pulse 
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(83) 


Suppose  the  vehicle  acceleration  is  described  by  a  half-sine 
pulse, 


m  =  -A.sm(y\0<,<r,, 

V 

where 


4  = 


7TV0 

4  Z>„ 


(72) 


(73) 


and  Tv  is  given  in  Eq.  (59).  The  velocity  and 
displacement  of  the  vehicle  are 

^(0  =  jy  +  ^^cos^,  (74) 

2  n  I, 


and 

*(0  =  jy/  +  j^^siny- 

[ 2  n  rv 

It  follows  from  Eqs.  (4),  (15),  and  (27)  that 

z  +  2£z  +  cozz  =  -^-(2^  o)+  co2t)  + 


(75) 


n 

Denote 

n 

C*V=T' 

1  V 

and 

a.  =  tan'<^). 

CO 

Then,  Eq.  (76)  becomes 
v, 


AJV  r_  J,  7T  f  2  •  "  •  n 

[2^rycos- —  +  cy  —sin — ] 


(76) 


T. 


n  T„ 


z  +  2£  ryz  +  <y2z  =  ®  +  <y20  + 


o.: 


l  +  (2£—  )2  sin(<yvf +  a2) 
I  co 


(77) 


(78) 


(79) 


By  solving  this  equation,  the  motion  of  the  occupant 
under  a  half-sine  pulse  is  found  to  be: 

2(t)  =  zc{t)  +  zpli). 
where 


zc(t)  =  Z sin  , 

and 

*,(o=^*+ 


(80) 

(81) 


A 

*7 


1  +  (2C— )2 

CD 


(82) 


sin(  &>v/  +  az  -  cp) 


[i-(— )2r+(^> 

I  CO  CO 

In  Eq.  (81), 


=  C>2  +CI  +  2C„Ca,ft> 
c  V  G-O®2 

and 


A  = tan"  ( 


,.cuVi7?-® 


). 


(84) 


2; 


where 


C,  — 

12  <y2 


ry. 


v\2 


1  +  (2C— ) 

co 


)!r+(— > 

I  ry  cy 


sin(a2-^) 


r  _vo  _A 

2 *  o  .^2 

2  ry, 


1  +  (2C— ) 

ry 


[l_(<yv)2]2+(2^)2 

ry  co 


cos  (a.  - 


In  Eqs.  (82)  and  (85), 

.  -t  2g(ocov 
^  =  tan  -2 — 2.. 

ry  -<y; 


(85) 

(86) 


The  occupant  responses  under  a  half-sine  pulse  are 
illustrated  in  Fig.  5,  where  co  =  32.78  and  £”  = 0.2 . 

Again,  the  optimal  Kinematics  of  the  occupant  is  not 
attained.  In  fact,  under  a  half-sine  crash  pulse,  it  is  also 
impossible  for  a  passive  restraint  system  to  provide  the 
required  protection  to  the  occupant  such  that  the  optimal 
kinematics  of  the  occupant  is  achieved  [8],  However, 
the  restraint  characteristics  can  still  be  optimized  to 
minimize  the  peak  occupant  deceleration  as  much  as 
possible.  Using  numerical  optimization  [8],  for  a  number 
of  levels  of  damping,  the  values  of  the  required 
rattlespace,  optimal  frequency,  and  the  corresponding 
minimum  peak  occupant  deceleration  are  listed  in  Table 
3. 


O  0  01  0  02  0.03  O.04  005  0.06  0.07  0.08  0  09  0  1 

Time  (s) 
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Figure  5.  The  motions  of  a  system  under  a  half-sine 
pulse 


Based  on  the  results  given  in  Table  3,  the  relationship 
between  the  peak  occupant  deceleration  and  the 
rattlespace  is  displayed  in  Fig.  3  by  the  red  curve.  The 
curve  indicates  that  the  peak  occupant  deceleration 
increases  with  the  increase  of  the  rattlespace,  because 
the  damping  is  controlled  at  each  level.  In  the  range  of 
[0.01  0.99],  larger  damping  results  in  smaller  required 
rattlespace  and  lower  peak  occupant  deceleration.  It  can 
be  seen  from  Fig.  3  among  the  three  crash  pulses 
considered,  the  optimal  crash  pulse  is  the  best,  followed 
by  the  constant-deceleration  pulse,  whereas  the  half¬ 
sine  pulse  is  the  worst,  in  terms  of  the  resulting  peak 
occupant  deceleration.  The  comparison  of  the  peak 
occupant  decelerations  between  the  optimal  crash  pulse 
and  the  half-sine  pulse  is  shown  in  Table  4  for  a  series 
of  specific  rattlespaces. 


Table  3.  Optimization  results  for  half-sine  pulse 


co  (rad/s) 

Am  (m/s2) 

S0(jn) 

0.01 

34.51 

414.76 

0.348 

0.10 

33.62 

365.16 

0.317 

0.20 

32.78 

331.12 

0.289 

0.30 

32.06 

310.20 

0.265 

0.40 

31.44 

297.02 

0.245 

0.50 

30.91 

288.46 

0.228 

0.60 

30.44 

282.75 

0.213 

0.70 

30.04 

278.87 

0.20 

0.80 

29.68 

276.18 

0.188 

0.90 

29.37 

274.24 

0.178 

0.99 

29.11 

272.96 

0.169 

Table  4.  Comparison  of  peak  occupant  decelerations 
_ between  optimal  pulse  and  half-sine  pulse 


S0  (m) 

(m/s2) 

(Optimal 

Pulse) 

Am  (m/s2) 

(Half-sine 

Pulse) 

Am  -Amo 

Amo 

x  1 00% 

0.169 

137.72 

272.96 

98.2 

0.178 

136.33 

274.24 

101.2 

0.188 

134.81 

276.18 

104.9 

0.20 

133.03 

278.87 

109.6 

0.213 

131.16 

282.75 

115.6 

0.228 

129.06 

288.46 

123.5 

0.245 

126.76 

297.02 

134.3 

0.265 

124.16 

310.20 

149.8 

0.289 

121.18 

331.12 

173.3 

0.317 

117.87 

365.16 

209.8 

0.348 

114.42 

414.76 

262.5 

CONCLUDING  REMARKS 


In  automobile  frontal  impact,  the  optimal  kinematics  of 
the  occupant  is  such  that  the  occupant  moves  at  a 
constant  deceleration.  For  given  restraint 
characteristics,  it  is  possible  for  the  occupant  to  attain 
optimal  kinematics  by  optimizing  the  vehicle  motion  or 
the  vehicle  crash  pulse.  For  linear  restraint 
characteristics,  the  optimal  vehicle  acceleration  should 
have  a  large  negative  impulse  at  the  initial  of  impact; 


after  that,  it  jumps  to  a  large  positive  value  and  then 
decays;  and  after  a  certain  period  of  time,  it  basically 
remains  constant  at  a  negative  value.  In  general,  the 
optimal  crash  pulse  depends  upon  the  characteristics  of 
particular  restraint  systems. 

Different  crash  pulses  have  different  influences  on  the 
occupant  motion  and  the  restraint  system  design.  For 
linear  restraint  characteristics,  the  peak  occupant 
deceleration  resulting  from  the  optimal  crash  pulse  is 
much  lower  than  those  produced  by  a  constant- 
deceleration  pulse  or  a  half-sine  pulse.  Under  the 
optimal  crash  pulse,  the  peak  deceleration  is 
independent  of  the  damping  of  restraint  systems. 
However,  certain  damping  in  the  systems  is  necessary, 
and  large  damping  can  reduce  the  amplitude  of  the 
initial  impulse.  Under  a  constant-deceleration  pulse  or  a 
half-sine  pulse,  the  damping  of  restraint  characteristics 
is  important  and  determines  the  required  rattlespace  and 
natural  frequency.  In  a  certain  range  of  the  damping,  the 
peak  occupant  deceleration  decreases  with  the  increase 
of  the  damping. 

It  is  recognized  that  the  use  of  a  lumped-parameter 
model  for  the  occupant-vehicle  system  and  the 
assumption  of  linear  characteristics  for  restraint  systems 
may  have  imposed  certain  limitations  on  the  analyses  in 
this  paper.  However,  the  results  and  conclusions  derived 
can  still  be  used  as  general  guidelines  for  the 
crashworthiness  design  of  the  vehicle  structure. 
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